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Growth of supermassive black holes and metallicity in quasars 
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Abstract. The strong correlation between the mass of the central supermassive black hole (SMBH) and the bulge in some 
galaxies and quasars implies that the formation of the black hole is somehow linked to the bulge. The measurement of metallicity 
by NV/CIV or NV/Hell in quasars allows to discuss a possible way of formation of the black hole. In this Letter we trace 
the metallicity along the possible routes in Rees' diagram in order to test the ways by which SMBHs can form. We derive a 
relation between the metallicity and the mass of the SMBH as Z oc Mbh based on the numerical simulation of the evolution 
of star clusters. It is in good agreement with the relation determined by the metallicity measured by NV/CIV or NV/Hell in 
Hamann & Ferland's sample. This lends observational support to the formation of SMBHs via routes R4 or R5, namely, the 
evolution of dense star clusters. 
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1. Introduction 

There is increasing evidence for the presence of super- 
massive black holes (SMBHs) in the centres of galax- 
ies and quasars from current observations (see the re- 
views of Rees 1984, Richstone et al 1998, Kormendy & 



Gebhardt 2001). The masses of SMBHs in galaxies strongly 
correlate with the masses of the corresponding bulges 
( |Magorrian et al 1998J), and even s tronger with the disper- 
sion velocity (Gebhardt et al. 2000). These facts reflect that 



the SMBHs grow with the bulges, albeit in a way which is 



yet unknown (Burkert & Silk 2001). The finding of the rela 



tion of Magorrian et al. (1998) may give us a clue to un- 
derstand how SMBHs are formed in galaxies. Laor (1998) 
originally realized that there is a similar relation between the 
masses of SMBHs and the bulge luminosities of their host 
galaxies in 14 bright quasars (for active galactic nuclei see 
Wandel 1999] ferrarese 'et al 200 1| ). In Laor's sample, SMBH 
masses are measured by the velocity dispersion of the H/3- 
emitting clouds, and the bulge luminosities are from Hubble 
Space Telescope (HST) observations. These relations are con- 
nected to the growth of the SMBHs. However, how exactly a 
black hole grows, remains open. 

Rees (1984) gave a map of all possible routes of SMBH 
formation. There are mainly two possible ways to form the 
central SMBH (see the routes of SMBH formation in Rees' 
diagram). The first is that SMBH forms directly from the pri- 
mordial gas cloud, this is further investigated by Loeb (1993), 
Loeb& Rasio (1994). The physical reason is that star formation 



would be quenched when the infalli ng primordial gas rea ches 
some critical central concentration (Loeb & Rasio 1994). On 



the other hand, a dense star cluster may be formed via star for- 
mation. This is supported by the observations of nearby galax- 
ies showing that there is a central star cluster with density of 
~ 10 6 — 10 8 A/q and a one dime nsional veloc ity dispersion a 
of typically ~ 100 - 400 Km s" 1 flLauer 1989|). The dense star 



cluster will inevitably form a SMBH, although the stars may 



evolve in somewhat different ways (Begelman & Rees 1978 
puncan & Shapiro 1983J |Quinlan & Shapiro 1987] )! 



It has been argued that three parameters, such as the mass 
and the spin of the black hole, and the accretion rate, may give 
a full description of a quasar (Blandford 1990). Recently, the 
metallicity, has received attention as the fourth parameter. It 
has become evident that quasars h ave very high metallicity , 
even at very high redshift (z > 4) ( Hamann & Ferland 1999 ). 
This means that there must have been very rapid and violent 
stellar evolution in th e early stages of these q uasars to produce 
the high metallicity (Haehnelt & Rees 1993). Metallicity, as a 
possible consequence of the formation and growth of SMBHs, 
may provide invaluable information of the history of quasars 
([Hamann & Ferland 1999). 
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We are motivated by the chemical evolution of quasars 
based on measurements of metallicity (Hamann & Ferland 
1992, 1993) to investigate this subject. The observation of 
emission line spectra of quasars can, in principle, provide the 
abundances in the broad line region (BLR), and would allow to 
deduce the details of the chemical evolution of the correspond- 
ing nucleus since the BLR is the innermost part of a quasar. In 
this Letter, we make an attempt to trace the metallicity of every 
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route in Rees' diagram, and naturally connect the growth of the 
SMBH to the observed metallicity in quasars. 

2. Observational metallicity in quasars 

High metallicity is common among quasars. As an indepen- 
dent parameter describing quasars, metallicity may be used to 
trace the evolutionary history of quasars. The metallicity in 
quasars may open a new way to understand the growth of a 
SMBH. It has been suggested by Hamann & Ferland (1993) 
that NV/CIV or NV/Hell are very good indicators of the 
metallicity in quasars. They subsequently measured the metal- 
licity in a large sample of quasars (Hamann & Ferland 1993, 
1999). There appears unambiguous evidence for a correlation 
between metallicity and luminosity. Although we do not fit this 
correlation via regression, it is found that 



or (Hi) a {vLv)C 



(i) 



where oc\ « a 2 ~ 0.5, and L v is the specific luminosity 
at frequency v. The reverberation technique provides a set of 
available data of quasars (Kaspi et al 2000), and shows a rela- 
tionship between the mass of a black hole and its luminosity 
Mbh oq (yhy) with « 0.5. From the above relations, we 
have 
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This relation is mainly derived observationally, and the only 
assumption is that the BLR clouds are virialized in the potential 
of the SMBH. It is generally believed that this approximation 
works in quasars (Netzer 1990). 

3. Metallicity in Rees' diagram 

Rees (1984) lists the possible ways to form SMBHs in the uni- 
verse. There are five possible routes. Route 1 (Rl) is collapse 
and/or accretion directly from an original gas cloud. Route 2 
(R2) is through post-Newtonian instability of a supermassive 
star. The supermassive star may be formed through two pos- 
sible ways: 1) direct collapse from a gas cloud, 2) from col- 
lisional disruption of stars in a dense star formed from the gas 
cloud. Route 3 (R3) is the merger of a black hole binary formed 
from the supermassive star induced by bar-mode instability. 
Route 4 (R4) means that a SMBH is born through spectacu- 
lar accretion in a cluster of neutron stars or stellar-mass black 
holes. Route 5 (R5) suggests that a SMBH is formed in a rela- 
tivistic cluster of compact objects as in R4, but the relativistic 
instability or gravitational radiation will induce the formation 
of SMBH. 

First we check route 1 (Rl). This route was studied by Loeb 
& Rasio (1994). In such a scenario, the SMBH is formed di- 
rectly due to a catastrophic collapse or it grows via gradual ac- 
cretion. As we show below, a supermassive black hole may be 
formed via a single catastrophic collapse, but without ejection 
of metals since the gravitational binding energy is too strong to 
eject matter. Thus there will be no relation between metallicity 
and the mass of the black hole. 



A supermassive star is first formed via two possible ways 
before collapsing to a SMBH. According to R2 a SMBH will 
be formed due to a post-Newtonian instability, R3 suggests that 
the SMBH is formed due to the merger of binary black holes. 
This process will lead to huge radiation of gravitational waves. 
It may be confirmed by the detection of gravitational radiation. 
The two routes (R2 and R3) employ the formation of a super- 
massive star. Here we will show that there is an upper limit for 
the mass of a SMBH created via R2. 

The structure and evolution of the supermassive star is in- 
sufficiently understood. We assume that the supermassive star 
is composed of hydrogen. The gravitational binding energy of 
an object with mass M and radius R is E g w GM 2 /R, where 
G is the gravitational constant. Its total nuclear energy is E n = 
eMc 2 with e being the conversion efficiency (e = 0.007 for the 
burning of hydrogen to helium) and c is the speed of light. If 
we adopt the mass-radius relation R = Ro (M/Mq) 11 for su- 
permassive stars, where Ro = 9.0 x 10 11 cm and q = 0.47 
(Collins 1989), we will have the upper limit of the mass of the 
supermassive object ejecting matter during the supernova ex- 
plosion 



^- = (2e) 1 /( 1 -«) (- 



2/(1-9) 



= 4.5 x 10 b 



(3) 



based on E g < E n , where vq — (GMq/Rq) ' . No mass 
ejection takes place during the formation of the black hole, al- 
beit the collapse of the supermassive star for M > M c . If a 
SMBH forms directly from the collapse of a primordial gas 
cloud, as the mass of SMBH is much higher than M c , the 
detectable metallicity may be low because (1) no heavy el- 
ements produced in the collapse have been ejected, and (2) 
the dominance of radiation pressure in supermassive stars and 
compact supermassive disks would prevent fragmentation and 
star formation in these systems (Loeb 1993, Loeb & Rasio 
1994). Although a strong wind may developed from the su- 
face of supermassive star, the first reason (1) is still sup- 
ported by the fact that the timescale of contraction of the su- 
permassive star is much shorter than that of the mass loss 
through the wind on the surface of supermassive star. The 
contraction timescale of a supermassive star approximates to 
^contr «3x 10 3 R U (M/M C ) 1/2 yr, where the t ypical dimen - 
sion of the supermassive star Ru = _R/10 14 cm ( Fricke 1973 ). 
This timescale is much shorter th an the mass loss timescal e 
t w = M c /M w 10 10 (M/Af c ) yrs ( |Bond, Arnet & Carr 1984| ). 
Therefore there is no significant ejection of metals during the 
rapid evolution of the supermassive star. 

It has been argued that the collision between a normal 
galaxy and a naked SMBH leads to the activity in a galaxy 
(Fukugita & Turner 1996), but in this case there is no correla- 
tion between metallicity and the mass of SMBH. On the other 
hand, if a SMBH is formed in a dense star cluster via coales- 
cence of stars or by a merger of stellar mass black holes, there 
is an approximate proportional relationship between the metal- 
licity and the mass of the SMBH according to the dynamical 
and stellar evolution of the cluster. We will derive a crude ap- 
proximation about this relation below. 
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Routes 4 and 5 (R4 and R5) suggest that the SMBH is 
formed in a dense star cluster. In R4 and R5, a cluster with 
stars heavier than lOOA/ is formed in advance. The high rate 
of supernova explosions leads to the formation of a cluster 
of neutron stars or steller-mass black holes. We assume that 
some of the ejected medium from the supernova explosions 
may form clouds in the broad line region. The two ways have 
a rough metallicity prediction as we show in the following. 
Following Duncan & Shapiro (1983), we assume that the clus- 
ter consists of N = 10 8 Ag stars with a dispersion velocity 
v = 350u35o(Kms _1 ), this star cluster thus has a core with a 
radius R c = 2.3Nsv^ pc. Out of the radius i? c , the initial den- 
sity has a power law profile n(r) ex r~ p , here p is the index. 
This power law distribution ma y work further into the klio par- 
sec scale bulge (Balcells 2001). The two-body relaxation time 



t; 



is t x = 1.5 X lO^iVgA^u^yr, where A i8 = ln(0.5iV)/18 
and the initial masses of all stars are taken as one solar mass 
m. The total mass and the collision time scale of star cluster 
are raN and 



t co ii = 6.8 x 10 10 iVJV 5 5 (l + 1.3^? 



:S5i) 



r 1 



yr. 



(4) 



respectively. The black hole grows due to the coalescence of 
stars via collision and tidal capture of stars because the col- 
lisional interactions are inelastic. Detailed numerical simula- 
tions show that the growth of a black hole can be divided into 
two phases for an isothermal sphere of solar type stars in the 
absence of a black hole (Duncan & Shapiro 1983): 



AWMeyr- 1 ) = 



9o, 

got 



for < t < t co n 
~ 7 , for t > t coU . 



(5) 



where go is a constant related to N& and W350, 7 = (3— 2p)/p w 
(0.9 ~ 1.0) is the growth index in the later phase weakly de- 
pending on the cluster parameters. For example, go « 1.7, 
7 = 1 and tcoix = 2.0 x 10 8 yr for the case of N 8 = 2.7, 
w 350 = 2.9 (Duncan & Shapiro 1983). Here we neglect the 
small peak of Mbh around the critical time t co \\. After the colli - 
sional phase the growth of the black hole is mainly via the cap- 
ture of stars from the cluster. It is thus expected that the main 
phase of metallicity enrichment is within the the time t co lh 

The Kelvin-Helmholtz time scale approximates to £kh ~ 
3.0 x 1Q 7 (m/ Mo)~ 2 yr, which represents the time scale to 
form a new star after coalescence (Quinlan & Shapiro 1987). 
Numerical results show that a massive star with about 50M Q 
can grow within a time scale t coa \ rs 10 6 yr due to colliding 
star coalescence if the dispersion velocity is less than its es- 
cape velocity (Colgate 1967). The maximum mass of a star due 
to coalescence via collision can be estimated from the fact that 
the coalescence process saturates when an ordinary IMq star 
cannot be captured by the coalescencing star since this star can 
pass straight through them. Supposing the mass-radius relation 
to be r*/R Q = (to*/M©) 3/4 , we have 



< 



v 

^csc 



M Q =53 



/«350\ f 

V 2.9/ 



Mr, 



(6) 



where v csc = (2GM q /Rq) 1 ^ 2 . This mass limit is sensitive to 
the dispersion velocity. Considering some other physical rules, 



the coalescencing stars may grow up to 50 ~ 100 solar mass 
( jBegelman & Rees 1978[ ). 

The time scale for star evolution approximates to t cvo — 
6.0 x 1O 6 (5OM /M)yr for stars larger than 12M , £ K h is 
much shorter than the coalescence and evolution time scales 
of massive stars. It appears that the solar mass stars in the 
cluster may form many massive stars following the collisions 
and coalescence within the collision time scale. Once a mas- 
sive star forms, the ejection of metal-rich ejecta will follow the 
post-main sequence nuclear burning, culminating most prob- 
ably in a Type II supernova. The heavy (Z > 6Z Q ) el- 
ement mass M z yield in ejecta varies with the progenitor 
mass, but is expected to be in the range of 4 ~ 40 M Q 
( Woosley & Weaver 1986 ). The rapid rotation of the progeni- 
tor can modify the metal yield by a modest amoun t as shown 
by the calculations ( Bodenheimer & Woosley 1983 ). Although 
it remains uncertain, the most likely remnants of massive pro- 
genitors are stellar-mass black holes with ~ 10 solar mass 
( Woosley & Weaver 1986 ). Considering that the life time of a 
massive star is much shorter than the collisional time scale, the 
total mass of metals produced by supernova explosion is about 
Mz = m z (dN/dt) t, then the metal abundance Z can be ob- 
tained by, 



V TO / \N dt J \ TO / Wcoll 



(7) 



Combining eqs (5) and (7) we obtain the following relation be- 
tween the mass of a growing black hole and the metallicity as 



I / 772 \ 

Z = — - — — Mbh, for < t < t colh 
g t C oi\ V m / 



(8) 



When t > toon, as the growth of the SMBH is dominated 
by the capture of stars, the above relationship will not hold. 
Detailed numerical calculations using a photoionization model 
show that the abundance Z is roughly proportional to the two 
ratios, 

.HeH J ' 



Z oc 



vcivy 



or Z ex 



(9) 



and 71 w 72 ss 1 measured from the right panel in Figure 6 
of Hamann & Ferland (1999). We would like to point out that 
this result is based on a series of broad emission line simula- 
tions that hold all other parameters fixed while the abundances 
are varied. We should keep this in mind. Thus the observable 
metallicity is given by 



^-\ ex M 1 /- 
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(10) 



This formula works within the collisional timescale t ca \\ which 
is determined by the cluster itself. We find that the theoretical 
prediction (10) is in good agreement with the observational re- 
lation (2). Therefore routes 4 and 5 (R4 and R5) may be the 
likely ways to form a SMBH from a primordial gas cloud. 

4. Conclusions and discussions 

In this paper we trace the metallicity along the routes in Rees' 
diagram in order to test the possible ways by which SMBHs 
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form in quasars. The argument that a single collapse of a su- 
permassive star forms a SMBH is not supported by the ob- 
served relation between metallicity and luminosity in quasars. 
The reason is that the mass is too high for the ejection of met- 
als since its self-gravitation energy is larger than the energy of 
nuclear reaction. We derive a formula of growth of a black hole 
with metallicity, which is in good agreement with the observed 
relation between metallicity and luminosity. Thus observations 
seem to support the formation of SMBHs from a dense star 
cluster via coalescence and evolution of stars. 

There is a modified route 1 (hereafter MR1) suggested by 
the referee. It states that if stars form during the collapse, they 
will enrich the gas with metals before the formation of the cen- 
tral SMBH. The gas that forms the SMBH might be processed 
by a generation of stars in the host galaxy before it gets fun- 
neled into the core to make the SMBH. Such a complex forma- 
tion of SMBH depends on many processes, fragmentation of 
cloud during the collapse, the formation of multi-generations 
of stars. Let us speculate that a primordial gas cloud with mass 
of several 10 8 M Q is gravitationally fragmentated into several 
hundreds of massive stars with mass less than 4.5xl0 6 M Q . 
These massive stars evolve undergoing hydrogen burning, then 
eject metals due to supernovae explosion, leaving a cluster with 
several hundreds of intermediate mass black holes (< 1O 6 M ). 
What's the observational consequence of such a cluster? This 
interesting route will be explored in the future. 

Recently, much atten tion has been given to th e recurrent ac- 
tivity model of quasars ( Haehnelt & Rees 1993|), w hich seems 
to be supported by observations ( Kuhn et al2001 ). The mass 
growth of a SMBH in the first generation of quasars may be 
mainly due to R4 or R5. The masses of SMBHs in recurrent 
quasars may not increase significantly in the subsequent gener- 
ations within the Eddingtontimescale A/BiiM^Edd sw 5.0 xlO 8 
years. This argument might be supported by the fact that the 
relation of the SMBH mass with the bulge mass in quasars 
(Laor 1998) is similar to the Magorrian et al. relation in galax- 
ies. We note that the slope of the relation found by Laor (1998) 
is slightly different from that of the Magorrian et al. relation. If 
the differences are real, we might then deduce the evolution of 
SMBH accretion after the first generation. This is worthy of a 
future investigation. 
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